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Abstract 

We study CP sensitive observables in neutralino production e + e~ — ► XiX^j an d 
the subsequent two-body decays of the neutralino Xi —> Xn% and of the Z boson 
Z — > U(qq). We identify the CP odd elements of the Z boson density matrix and 
propose CP sensitive triple-product asymmetries. We calculate these observables 
and the cross sections in the Minimal Supersymmetric Standard Model with complex 
parameters \x and M\ for an e + e~ linear collider with y/s = 800 GeV and longitu- 
dinally polarized beams. We show that the asymmetries can reach 3% for Z — > it 
and 18% for Z — > qq and discuss the feasibility of measuring these asymmetries. 
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1 Introduction 



In the Minimal Supersymmetric Standard Model (MSSM) [1] several supersymmetric 
(SUSY) parameters can be complex. In the neutralino sector of the MSSM these are 
the U(l) gaugino mass parameter Mi and the Higgsino mass parameter /i. (The SU(2) 
gaugino mass parameter M2 can be made real by redefining the fields.) The physical 
phases ipM 1 and (p^ of Mi and //, respectively, imply CP odd observables which can in 
principle be large, because they are already present at tree level. It has been shown that 
in the production of two different neutralinos e + e~ — > Xi'X.'j the CP violating phases cause 
a non- vanishing neutralino polarization perpendicular to the production plane [2, 3, 4, 5], 
which leads to CP odd triple-product asymmetries [6] of the neutralino decay products 
[4, 5, 7, 8, 9]. 

In this work we study CP violation in neutralino production 

e + + e- -> + i,j = l,...,4, (1) 

with the subsequent two-body decay of one neutralino into the Z boson (for recent studies 
see [9, 10]) 



r,0 



X 



Xn + Z; n<t, (2) 



and the decay of the Z boson 

Z->f + f, f = £,q, £ = e,fi,T, q = c,b. (3) 

In case of CP violation the non-vanishing phases ipMi and tp^ lead to CP sensitive elements 
of the Z boson density matrix, which we will discuss in detail. Moreover, these CP 
sensitive elements cause CP odd asymmetries Af in the decay distribution of the decay 
fermions [4]: 

A f = a(Tf >0) ~ a(Tf < 0) (4) 
f ' a(T f > 0) + a(T f < 0) ' 1 ; 

with a the cross section and the triple product 

Tf = Pe- ■ (Pf x Pf). (5) 

Due to the correlations between the polarization and the Z boson polarization, there 
are CP odd contributions to the Z boson density matrix and to the asymmetries from the 
production (1) and from the decay process (2). 

The triple product Tf, Eq. (5) changes sign under time reversal and is thus T odd. 
Due to CPT invariance, the corresponding T odd asymmetries Af are also CP odd if final 
state interactions are neglected. The final state interactions would also contribute to Af. 
However, they only arise at loop level and are neglected in the present work. 

In Section 2 we give our definitions and the formalism used and define the Z boson 
density matrix. In Section 3 we discuss some general properties of the asymmetries. We 
present numerical results in Section 4. Section 5 gives a summary and conclusions. 
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2 Definitions and formalism 



We give the analytic formulae for the differential cross section of neutralino production 

e + + e" -> XiiPxi'^ + xSiPxj^j), (6) 

with longitudinally polarized beams and the subsequent decay chain of one of the neu- 
tralinos 

X°i ^ X°n(p Xn An) + Z(p Z ,\ k ), (7) 

Z -> f(p f ,X f )+f(p f ,X f ). (8) 

In Eq. (6) and Eq. (7), (8), p and A denote momentum and helicity, respectively. For 
a schematic picture of the neutralino production and decay process see Fig. 1. In the 
following we will derive the Z boson spin-density matrix and relate it to the CP asymmetry 
A f in Eq. (4). 




Figure 1: Schematic picture of the neutralino production and decay process. 



2.1 Lagrangian and helicity amplitudes 

The interaction Lagrangians relevant for our study are (in our notation and conventions 
we follow closely [1, 9]): 



C 



L 
C 



\Z»X^[OjP L + 0*P R }xl i,j 
gf^P R xTe L + 9f*eP L X-e R + h.c, 
Zj^[L f P L + R f P R ]f, 



,4, 



'Z°ff 

with Pl,r = |(1 T 75)- I n the neutralino basis 7, Z, H®, H® the couplings are: 



O 
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ij 



2 cos 9w 



(N l3 N* 3 - N t4 N* A ) cos 1(5 + (iV l3 jV* 4 + N U N* 3 ) sin 2/3 



(9) 
(10) 

(11) 
(12) 



3 



o 



'R 



-0" L * 



9 



-V2 



1 



cos 9 W 

2 



(T 3f - q f sm 2 9 W ), R f 



cos 9 



w 



-qf sin 2 9 W , 



(T 3i - q e sin 9 w )N i2 + q e sin 9 w N n 



cos9 w 

fg = -V2q e sin 9 W [tan 9 W N* 2 - iV* 



(13) 

(14) 
(15) 



with g the weak coupling constant (g = ej sin^, e > 0), qf and T 3 j the charge and the 
isospin of the fermion, and tan/5 = i^/fi the ratio of the vacuum expectation values of 
the two neutral Higgs fields. is the complex unitary 4x4 matrix which diagonalizes 
the neutral gaugino-Higgsino mass matrix Y a/3 , N* a Y a/3 N^ k = m x o5ik, with m x q > 0. Note 



that our definitions of 0^ ,R and Lf,Rf differ from those given in [1, 3] by a factor of 

gj COs9 W . 

The helicity amplitudes Tp Xj for the production process are given in [3]. Those for 
the two-body decays, Eq. (7) and Eq. (8), are 



1 D 1 ,X i 



(16) 



and 

rp X f X f 
£>2,Afc 



: u( Pf ,\ f )r[L f P L + R f P R Hp f -,X f )e^. (17) 

The polarization vectors e^, = 0, ±1, are given in Appendix A. The amplitude for the 
whole process (6), (7), (8) is 

T = A(x?)A(Z) £ T^T^^l, 



A^A^ 



(18) 

m 2 . + im x T Xi ] and the Z boson prop- 
agator A(Z) = i/\pz ~ m z + imz^z] (the mass and width are denoted by m and T, 
respectively). For these propagators we use the narrow width approximation. 



with the neutralino propagator A(x°) = 



2.2 Cross section and Z boson density matrix 

For the calculation of the cross section for the combined process of neutralino production 
(6) and the subsequent two-body decays (7), (8) of Xi we use the same spin-density matrix 
formalism as in [3, 11]. The (unnormalized) spin-density matrix of the Z boson 

p P (Z)W = |A(x°)| 2 E Pp(x!) AlA 'P fl i(x?)^, (19) 

A;,A^ 

is composed of the spin-density production matrix 
and the decay matrix 

PnM) X $ = Y, T kx T DuK- ( 21 ) 

An 
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With the decay matrix for the Z decay 

pD2(Z) X ' k X k = J2 T dIX T D 2 X ( 22 ) 



the amplitude squared for the complete process e + e — > XiX^] Xi ~^ XnZ', Z — > / / can 
now be written 

|T| 2 = \A(Z)f ^ Pp(Z) XkK PD2(Z) KXk . (23) 

The differential cross section in the laboratory system is then given by 

dcr= ^|T| 2 rfLips(s,p Xj ,p Xn ,p / ,p / -), (24) 

where dLips(s,p Xj ,p Xn ,pf,pf) is the Lorentz invariant phase space element defined in 
Eq. (B.l) of Appendix B. More details concerning kinematics and phase space can be 
found in Appendices A and B. 

For the polarization of the decaying neutralino Xi with momentum p Xi we introduce 
three space like spin vectors s x . (a = 1,2,3), which together with p^./m Xi form an or- 
thonormal set with s x . ■ s b x . = —5 ab , s£. • p Xi = 0, then the (unnormalized) neutralino 
density matrix can be expanded in terms of the Pauli matrices: 

p P (x°) AiA * = 2(5 A ^P + ^ A ,£ P ), (25) 

where we sum over a. With our choice of the spin vectors s£\, given in Appendix A, -g- 
is the longitudinal polarization of neutralino xl, ~p" is the transverse polarization in the 

E 2 

production plane and is the polarization perpendicular to the production plane. The 
analytical formulae for P and E P are given in [3]. To describe the polarization states of 
the Z boson, we introduce a set of spin vectors t° z (c = 1,2,3) and choose polarization 
vectors s^ k (\k = 0, ±1), given in Appendix A. Then we obtain for the decay matrices 

PnMf$ = (^^r + 4A,^r)^^ (26) 
and 

PD2(Z) KXk = D^e^\ (27) 
with [9]: 

D? = 2[2p^ Xi -(p^+p x ^)-|(m Xj +m Xn -m|)^]|0:f| 2 

-2g^m Xi m Xn [{ReO'^f - (ImOj) 2 }, (28) 
£ a D f = 2i{-m Xi 6^s x Jp Xi p - pz,)\0'^\ 2 + 2m Xn (s a x y xi - s^){ImO^){ReO^) 
-m Xn e^ s«p XiP [{ReO^Y - (ImOj) 2 ]}; (e 0123 = 1), (29) 
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and 

= 2(-2p^- + p^- + p^-§m|^ (30) 

Due to the Majorana properties of the neutralinos, is symmetric and £^{ tiy is anti- 
symmetric under interchange of p and v. In Eq. (26) and Eq. (27) we use the expansion 
[12]: 

e x *et' k * = |^^V-^e^ CT p^(J c ) A '^-^| M 4,(^) A ^ fc , (^0123 = 1), (31) 
summed over c, d. Here, J c are the 3x3 spin 1 matrices with [J c , J d ] = ie C( i e J e and 

jcd = Jcjd + J djc_4 S cd^ ^ 

with J 11 + J 22 + J 33 = 0, are the components of a symmetric, traceless tensor, given in 
Appendix C, and 

T VZjxVZv (V3 s 

L iiv — 9p,v -r 2 [o6) 
m z 

guarantees the completeness relation of the polarization vectors 

£^ = + T (34) 

x k m z 

The second term in Eq. (31) describes the vector polarization and the third term describes 
the tensor polarization of the spin 1 Z boson. The decay matrices can be expanded in 
terms of the spin matrices J c and J cd . The first term of the decay matrix p^, Eq. (26), 
which is independent of the neutralino polarization, then gives 

D^e Xk *Su k = Di<J AfcA * + c D 1 (J c ) XkX 'k + cd D 1 (J cd ) x " x 'i<, (35) 
summed over c, d, with 

.2 12 _2 , ^{Pxr-Pzf- 



D 1 = 



m v - ^m v . - m z + T 

L\2 lT^n,"L\2 



L|2 



+2m Xi m Xn [(ReO n tY - {ImO n ^)% (36) 



cei 



L\2 



Di = -[m-V x M-V Xl ) + \{rn 2 Xi +m 2 Xn -m 2 z )5^ 

-5 cd m Xi m Xn [(ReO:ty - (ImOj) 2 }, (37) 

and C D\ = due to the Majorana character of the neutralinos. As a consequence of the 
completeness relation, Eq. (34), the diagonal coefficients are linearly dependent 

U D 1 + 22 D 1 + 33 D 1 = -§£>!. (38) 

For large three momentum p Xi , the Z boson will mainly be emitted into the forward 
direction with respect to p Xi , i.e. p Xi ~ pz, with p = p/\p\, so that (t z 2 ■ p Xi ) ~ 
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in Eq. (37). Therefore, for high energies 11 D\ 22 Di, and the contributions for the 
non-diagonal coefficients cd Di(c=^ d) will be small. 

For the second term of pm, Eq. (26), which depends on the polarization of the decaying 
neutralino, we obtain 



^Dl V tv 



summed over c, d, with 



(39) 



C £di = ^{[\0:f\ 2 m Xi + [(ReO:f) 2 -(ImOj) 2 }m Xri 
m z >■ 



x 



X • PzWz -Pxi)- (4 ■ tz)(Pz ■ Pxd + \O n f\ 2 m Xi m 2 z (s Xi ■ t%) 



-2(7mOl L )(it:eOl L )m x „e^ CTS K i p^}, 



(40) 



and EJ )1 = E^ )1 = due to the Majorana character of the neutralinos. A similar 
expansion for the Z decay matrix, Eq. (27), results in 



PD2(Z) KXk = D 2 5 x '* Xk + c D 2 (J c ) x ' kXk + cd D 2 (J cd ) x ^ Xk , 
where we sum over c, d, with 



D 2 = l(Rj + L})m 2 z , 
C D 2 = 2(R 2 f -L 2 f )m z (t z -p f ), 
<*D 2 = {R) + L))W z . Vj ){t d z . Vj )-\m 2 z ^ 



(41) 

(42) 
(43) 
(44) 



As a consequence of the completeness relation, Eq. (34), the diagonal coefficients are 
linearly dependent 



ii 



D 2 + 22 D 2 + 4i D 2 = 



33 



3 n 



(45) 



For large three-momentum p z , the fermion / will mainly be emitted into the forward 
direction with respect to p z , i.e. p z ~ pj, so that {t l z 2 • pj) ps in Eq. (44). Therefore, 
for high energies ll D 2 ps 22 -D2, and the contributions for the non-diagonal coefficients 
cd D 2 {c=f d) will be small. 

Inserting the density matrices (25) and (26) into Eq. (19) leads to: 



p P (Z) XkX 'k = 4 |A(x-)| 2 [PD 1 5 XkX * + E^ C E^ (J C ) A ^ + P cd D 1 {J cd ) XkX 'k 
summed over a, c, d. Inserting then (46) and (27) into Eq. (23) leads to: 



(46) 



\T\ 2 = 4 |A(x-)| | A(Z) | 2 [3P£>i£>2 + 2E£ c T, a D1 C D 2 + AP( cd D cd D 2 - § CC D 1 dd D 2 ) , (47) 

summed over a, c, <i, which is the decomposition of the amplitude squared in its scalar 
(first term), vector (second term) and tensor part (third term). 
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2.3 Z boson density matrix 

The polarization of the Z boson, produced in the neutralino decay (7), is given by its 
3x3 density matrix < p(Z) > with Tr{< p(Z) >} = 1. We obtain < p(Z) > in the 
laboratory system by integrating Eq. (46) over the Lorentz invariant phase space element 
dLips(s,p X] ,p Xn ,p z ) = dLips(s,p Xi ,p Xj .) ds Xi £ ± dLips(s Xi ,p Xn ,pf), see Eq. (B.l), 
and normalizing by the trace: 



summed over c, d. The vector and tensor coefficients V c and are given by: 



(48) 



f|Afe°)| 2 ^ P c S a D1 rfLips _ /|A(x°)| 2 P cd A ^Lips 

C 3/|A(x°)| 2 P^rfLips ' ±cd ±dc 3/|A(x°)| 2 P^rfLips' ^ 

with sum over a. The tensor coefficients T\ 2 and T23 vanish due to phase-space integration. 
The density matrix in the circular basis, see Eq. (A. 11), is given by 

(50) 
(51) 
(52) 

(53) 
(54) 

where we have used Tn + T 22 + T 33 = — \ and T i2 = T23 = 0. 





> 


= I-V3 + T33, 


< P(Z) 00 


> 


= -2T33, 


< p(z)-° 


> 


= j=(V 1 + tV 2 )-V2T 13 


< p{zy + 


> 


= T n , 


< p(zr 


> 


= j= 2 (V 1 + iV 2 ) + V2T 13 



3 T odd asymmetry 

From Eq. (47) one obtains for the asymmetry, Eq. (4): 

JSignpT/llTfdLips = /|A(x°)| 2 |A(Z)| 2 Sign[T / ]2Sg ) ^ 1 ^Lips 
f J|T| 2 rfLips J|A(x°)| 2 |A(Z)| 2 3PZnr; 2 dLips ' 1 ' 

summed over a, c. In the numerator only the vector part of |T| 2 remains because only 
the vector part contains the triple product 1 7} = p e - ■ (pf x pj). In the denominator 
the vector and tensor parts of |T| 2 vanish, because for complete phase space integrations 
the spin correlations are eliminated. Due to the correlations between the and the Z 
boson polarization, Sp c SJ )1 , there are CP odd contributions to the asymmetry Af which 
stem from the neutralino production process, see Eq. (6), and/or from the neutralino 
decay process, see Eq. (7). The contribution from the production is given by the term 
with a = 2 in Eq. (55) and it is proportional to Sp, Eq. (25), which is the transverse 

^Note that if one would replace the triple product 7} by 7} = p e - ■ (p Xi x pz), and would calculate 
the corresponding asymmetry, where the Z boson polarization is summed, all spin correlations and thus 
this asymmetry would vanish identically because of the Majorana properties of the neutralinos. 
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polarization of the neutralino perpendicular to the production plane. For e + e~ — > xlXi 
we have T? P = 0. The contributions from the decay, which are the terms with a = 1, 3 in 
Eq. (55), are proportional to 

c X a D1 c D 2 => -8m x J/mO:f)(it>eO:f)^ (56) 

see last term of Eq. (40), which contains the e-tensor. Thus Af can be enhanced (reduced) 
if the contributions from production and decay have the same (opposite) sign. Note that 
the contributions from the decay would vanish for a two-body decay of the neutralino 
into a scalar particle. In this case the remaining contributions from the production are 
multiplied by a decay factor oc (|P| 2 - \L\ 2 ) [7], and thus Af oc (\R\ 2 - \L\ 2 )/(\R\ 2 + \L\ 2 ), 
where R and L are the right and left couplings of the scalar particle to the neutralino. 

For the measurement of Af the charges and the flavors of / and / have to be distin- 
guished. For / = e, p this will be possible on an event by event basis. For / = r it will be 
possible after taking into account corrections due to the reconstruction of the r momen- 
tum. For / = q the distinction of the quark flavors should be possible by flavor tagging 
in the case q = b,c [13]. However, in this case the quark charges will be distinguished 
statistically for a given event sample only [14]. Note that A q is always larger than Ae, 
due to the dependence of Af on the Z-f-f couplings [4, 9]: 

^ -V> = 5 + H ffi" I ffi ■* * 6.3 (4.5) x At, (57) 

which follows from Eqs. (42), (43) and (55). 

The relative statistical error of Af is given by 5Af = AAf/\Af \ = Sf/(\Af\y/~N) [7], 
with Sf standard deviations and N — C ■ a t the number of events with C the integrated 
luminosity and the cross section a t = a(e + e~ — > x BR(x° — > Zx° n ) x BR(Z — > //). 

Taking 5Af = 1 it follows Sf = \Af\y/~N. Note that Sf is larger for f = b,c than for 
/ = £ — e,/i,T with S b ~ 7.7 x S e and S c ~ 4.9 x S e , which follows from Eq. (57) and 
from BR(Z -> 66) ~ 1.5 x BR(Z -> ££), BR(Z -> cc) ~ 1.2 x BR(Z -> il). 

4 Numerical results 

We present numerical results for the Z density matrix < p(Z) >, Eq. (48), the asymmetry 
Ae(£ = e,fx,r), Eq. (4), and the cross section a t = a(e + e~ — > XiX°j) x BR(x° — >■ X?^) x 
BR(Z — > ££). For the branching ratio Z — > summed over i = e,/j,,r, we take the 
experimental value BR(Z — > ££) = 0.1 [15]. The values for ^4^ ,, may be obtained from 
Eq. (57). We choose a center of mass energy of ^/s = 800 GeV and longitudinally polarized 
beams with beam polarizations (P e -,P e +) = (±0.8, =f0.6). We study the dependence of 
< p(Z) >, At and (T t on the MSSM parameters p = \p\e iLp » and M x = \M x \e iLpM i. 
For all scenarios we keep tan/5 = 10. In order to reduce the number of parameters, we 
assume the relation \M\\ = 5/3M2tan 2 #jy and use the renormalization group equations 
[16] for the selectron and smuon masses, m| = ijiq + 0.23Mf — m 2 z cos 2/3 sin 2 9w, to- = 
m 2 + 0.79M 2 ±m|cos2/5(-l/2±sin 2 ^), taking m = 300 GeV. 
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Our numerical results presented below are obtained at tree level. One-loop corrections 
to e + e~ — > x°iX°j nave been given in [17] for real MSSM parameters. They are of the order 
of a few percent and may reach values up to 10%. As the bulk of the one-loop corrections 
are presumably CP-even, we expect that they will not significantly change our tree-level 
result for Af. For an appropriate analysis of the one-loop corrections to Af it would be 
necessary to adopt the formulae of [17] to the case of complex MSSM parameters, which 
is beyond the scope of the present paper. 

The experimental upper limits on the electric dipole moments (EDMs) of electron 
and neutron may restrict the phases y? M and <^m 1 - These restrictions are very model 
dependent. They are less severe when cancellations between the contributions of different 
SUSY phases occur. For example, in the constrained MSSM the phase y?^ is restricted 
to \<Pn\ < 0.17T, whereas the phase y^ is not restricted, but correlated with y?^ [18]. In 
most of our numerical examples below we have chosen ip>M 1 = ±tt/2, y? M = 0, which agrees 
with the constraints from the electron and neutron EDMs. In order to show the full 
phase dependences of the CP asymmetry Af, in one example we study its y?^ behavior 
in the whole y?^ range, relaxing in this case the restrictions from the EDMs. However, as 
shown in [19], if also lepton flavor violating terms are included, the EDM constraints on 
<ffj, disappear. 

For the calculation of the neutralino widths T Xi and the branching ratios BR(x? — > 
X\Z) we neglect three-body decays and include the following two-body decays, if kine- 
matically allowed, 

Xi -> e RyL e, j2 RiL ii, f m r, v t v h x° n Z, Xm W± , X° n H u ^ = e,/i,r, to = 1,2, n<i (58) 

with being the lightest neutral Higgs boson. The Higgs parameter is chosen uia = 
1000 GeV and thus the decays Xi ~^ xt,H T into the the charged Higgs bosons, and the 
decays Xi ~^ Xn H23 i n t° the heavy neutral Higgs bosons are forbidden in our scenarios. 
In the stau sector, we fix the trilinear scalar coupling parameter A T = 250 GeV. 

4.1 Production of X1X2 

In Fig. 2a we show the cross section for X1X2 production in the \/j,\-M 2 plane for = 
and <p>M 1 — 0.57r. For > 250 GeV the left selectron exchange dominates due to 
the larger X2 ~~ &l coupling, so that the choice of polarization (P e -,P e +) = (—0.8,0.6) 
enhances the cross section, which reaches values of more than 110 fb. The branching 
ratio BR(x° — ► -^Xi) is shown in Fig. 2b. The branching ratio can even be 100% and 
decreases with increasing and M 2 , when the two-body decays into sleptons and/or 
into the lightest neutral Higgs boson are kinematically allowed. The cross section a t = 
a(e + e~ -> xXxl) x BR (X2 ~? z Xi) x B R(^ -»• ®) is shown in Fig. 2c. Due to the small 
branching ratio BR(Z — > it) = 0.1, a t does not exceed 7 fb. Fig. 2d shows the \/j,\-M 2 
dependence of the asymmetry Ae for = 0.57T and ip^ = 0. The asymmetry \Ae\ can 
reach a value of 1.6%. On the contour in Fig. 2d, the (positive) contributions from the 
production cancel the (negative) contributions from the decay. We also studied the y^ 



10 



dependence of At- In the \p\-M 2 plane for ipM 1 = and ip^ = 0.5n we found \Ae\ < 0.1 



In Fig. 3 we show the vfy-VMi dependence of At for = 400 GeV and M 2 = 250 GeV. 
The value of Ae depends stronger on (/^i than on tp^. It is remarkable that the maximal 
phases of (pM 1 , i Pn = ±tt/2 do not lead to the highest values of Ae ~ ±1.4%, which are 
reached for (ifiMi, w) ~ (±0.37r, 0). The reason for this is that the spin-correlation terms 
Ep c E^ 1 c D 2 in the numerator of Af, Eq. (55), are products of CP odd and CP even 
factors. The CP odd (CP even) factors have a sine-like (cosine-like) phase dependence. 
Therefore, the maximum of the CP asymmetry Af is shifted from (fM^f^ = ±7r/2 to a 
smaller or larger value. 

In the (Pfx-ifiMi region shown in Fig. 3 also the cross section a t = <r(e + e~ — > X1X2) x 
BR(x° 2 -> Zxl) x BR(Z -> if) with BR( X ° 2 -> = 1 and BR(Z -> it) = 0.1, is 

rather insensitive to and ranges between 7 fb (</?Mi = 0) and 14 fb (</?Mi = ±vr). 
For the leptonic decay of the Z, the standard deviations are given by St = \Ae\y/C ■ &t, 
and for the hadronic decays by S^ c ) = 7.7(4.9)Sg, see Section 3. For C = 500 fb -1 and 
(<fM 1 , ( P fJ ) = (±0.37r,0) in Fig. 3 we find Sb( c ) = 8(5) and thus Ab( c ) could be measured. 
However note that we have Se < 1 in this scenario and thus Ae cannot be measured at the 
68% confidence level (Se = 1)- In Fig. 4 we show the lpm x dependence of the vector (Vi) 
and tensor (To) elements of the Z density matrix < p(Z) >. The elements Tn, T22 and 
V\ have a CP even behavior. The element V2 is CP odd and is not only zero at y?Mi = 
and (fMi = 7r, but also at v?Mi ~ (1 ± 0.2)7r, which is due to the destructive interference 
of the contributions from CP violation in production and decay. The interference of the 
contributions from the CP even effects in production and decay cause the two maxima 
of V\. As discussed in Section 2.2, the tensor elements T u and T 22 are almost equal. 
Compared to V\ and V 2 , they have the same order of magnitude but their dependence 
on (fM 1 is rather weak. Furthermore, the other elements are small, i.e. Ti3,V3 < 10~ 6 
and thus the density matrix < p(Z) > assumes a symmetric shape. In the CP conserving 
case, e.g. for <p Ml = ^ = 0, M 2 = 250 GeV, \p\ = 400 GeV, tan/3 = 10, m = 300 GeV, 
v/i = 800 GeV and (P e -,P e +) = (-0.8,0.6) it reads: 



< P(Z) >= 



( 0.329 0.049 0.0003 

0.049 0.343 0.049 | . (59) 
V 0.0003 0.049 0.329 



In the CP violating case, e.g. for ifM 1 = 0.5n and the other parameters as above, < p(Z) > 
has imaginary parts due to a non- vanishing V^: 



< P(Z) >- 



I 0.324 0.107 + 0.037i 0.0003 

0.107 - 0.037i 0.352 0.107 + 0.037i |. (60) 
V 0.0003 0.107 -0.037i 0.324 



Imaginary parts of < p(Z) > are thus an indication of CP violation. Note that also the 
diagonal elements, being CP even quantities, are changed for ipu 1 =f and y? M y 0. This 
fact has been exploited in [10] as a possibility to determine the CP violating phases. 
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4.2 Production of X2 x\ 

In Fig. 5a we show the cross section a t = a(e + e~ — > X2X2) x BR(X2 — > ^Xi) x BR(Z — > 
in the \/j,\-M 2 plane for y? M = and pu x — 0.57T. The production cross section cr(e + e~ — > 
X2X2)' which is not shown, is enhanced by the choice (P e -,P e +) = (—0.8, 0.6) and reaches 
values up to 130 fb. The branching ratio BR^ — ► ^X?), shown in Fig. 2b, can be 100%. 
However, due to the small branching ratio BR(Z — > ££) = 0.1, the cross section shown in 
Fig. 5a does not exceed 13 fb. 

If two equal neutralinos are produced, the CP sensitive transverse polarization of the 
neutralinos perpendicular to the production plane vanishes, Ep = in Eq. (55). However, 
the asymmetry Af need not vanish, because there are CP sensitive contributions from the 
neutralino decay process, terms with a = 1, 3 in Eq. (56). In Fig. 5b we show the and 
M 2 dependence of the asymmetry Ae, which reaches more than 3% for </?Mi = 0.57T and 
(p^ — 0. Along the zero contour in Fig. 5b the contribution to At which is proportional 
to Ep, see Eq. 55, cancels that which is proportional to Ep. As the largest values of 
Ag > 0.2% and A q > 1% lie in a region of the \n\-M 2 plane where a t < 0.3 fb, it will 
be difficult to measure Af in a statistically significant way. We also studied the p^ 
dependence of Ae- In the \fj,\-M 2 plane for p>M 1 = and ip^ = 0.5n we found \Ae\ < 0.5%, 
and thus the influence of ip^ is also small. 

In Fig. 6 we show the pu x dependence of the vector (Vi) and tensor (T u ) elements of 
the Z density matrix < p(Z) >. Because there are only CP sensitive contributions from 
the neutralino decay process, V2 is only zero at pu x = 0,n and V\ has one maximum at 
Pm 1 = compared to the elements shown in Fig. 4. In addition, in Fig. 6 the vector 
elements V\ and V2 are much smaller than the tensor elements Tu ~ T22, compared to 
Fig. 4. The smallness of the vector element V2 accounts for the smallness of the asymmetry 
\Ae\ < 0.05%. Furthermore, the other elements are small, i.e. T 13 < 10~ 6 and V 3 = 0. 

4.3 Production of X\ X3 

In Fig. 7a we show the cross section a t = a(e + e~ — > X1X3) x BR(X3 —> ZXi) x BR(Z — > ££) 
in the \/i\-M 2 plane for </? M = and pu Y = 0.57T. The production cross section <r(e + e~ — > 
XiX-j), which is not shown, is enhanced by the choice (P e -, P e +) = (0.8, —0.6) and reaches 
up to 50 fb. The branching ratio BR(xg — > Zxi), which is not shown, can be 1. However, 
due to the small branching ratio BR(Z — > £€) = 0.1, the cross section shown in Fig. 7a 
does not exceed 5 fb. In Fig. 7b we show the \/j,\-M 2 dependence of the asymmetry Ae- 
The asymmetry \Ae\ reaches 1.3% at its maximum, however in a region, where a t < 0.3 
fb, the asymmetry Ae thus cannot be measured. We also studied the p^ dependence of 
Ae- In the \p\-M 2 plane for pu 1 = and p^ = 0.5n we found \Ae\ < 0.7%. 

4.4 Production of x\ X3 

For the process e + e~ — > X2X3 we discuss the decay X3 ~ * Zx± of the heavier neutralino 
which has a larger kinematically allowed region than that of x\ ~~ * Zx±- Similar to X1X3 
production and decay, the cross section a(e + e _ — > X2X3) reaches values up to 50 fb for a 
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beam polarization of (P e -,P e +) = (0.8, —0.6). The cross section for the complete process 
a t = a(e+e- -> X2X3) x BR(x° -> ^X?) x BR(Z -> attains values up to 5 fb in the 
\/j,\-M 2 plane, see Fig. 8a. 

The asymmetry Ae, Fig. 8b, is somewhat larger than the asymmetry for x® \\ produc- 
tion and decay, and reaches at its maximum 2%. Although in the respective region the 
cross section is also a bit larger, o t < 4 fb, it will be difficult to measure Ae- For example 
taking \fj,\ = 380 GeV, M 2 = 560 GeV and (</? Ml , </? M ) = (0.57T, 0), we have 1, for 

C = 500 fb -1 . However for the hadronic decays of the Z we have S b ( c ) 8(5) and thus 
Ab( c ) could be measured for X1X3 production. Concerning the ip^ dependence of Ae we 
found that \Ae\ < 1% in regions of the \y1\-M2 plane where at < 0.5 fb, and \Ae\ ^ 0.4% 
in regions where a t < 5 fb, for example for ip^ = 0.57T and y?Mi = 0. 

5 Summary and conclusions 

We have proposed and analyzed CP sensitive observables in neutralino production e + e~ — > 
X^,X°j an d the subsequent two-body decay of one neutralino into the Z boson ~^ Xn^, 
followed by the decay Z — > 11 for I = e, //, r, or Z — > for q — c,b. The CP sensitive 
observables are defined by the vector component V2 of the Z boson density matrix and the 
CP asymmetry Ae( q ), which involves the triple product %( q ) = p e - • (pe(q) x Pe(q))- The tree 
level contributions to these observables are due to correlations of the neutralino xi s Pi n 
and the Z boson spin. In a numerical study of the MSSM parameter space with complex 
Mi and \x for X1X21 X2X2) X1X3 an d X2X3 production, we have shown that the asymmetry 
Ae can go up to 3%. For the hadronic decays of the Z boson, larger asymmetries are 
obtained with A c (b) — 6.3(4.5) x Ae- By analyzing their statistical errors, we found 
that the asymmetries A c (b) could be accessible in future electron positron linear collider 
experiments in the 500-800 GeV range with high luminosity and longitudinally polarized 
beams. 
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Appendix 



A Coordinate frame and spin vectors 

We choose a coordinate frame in the laboratory system such that the momentum of 
neutralino Xj points in the z-direction (in our definitions we follow closely [3]). The 
scattering angle is l{p e - , p Xj ) and the azimuth can be chosen zero. The momenta are 
given by: 

p e - — E b (l, — sin 0,0, cos0), p e + = E b (l, sin0, 0, — cos0), (A.l) 
= 0,0, -<z), p Xj = (E Xj ,0,0,q), (A.2) 
with the beam energy E b = \/s/2 and 

s + m 2 — m 2 . s + ml —ml A^(s, ml , ml.) 

E x . = E x . = ^, q= . Xj , (A.3) 

where m Xi ,m x . are the masses of the neutralinos and X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + xz + 
yz). We choose the three spin vectors s x f (a = 1, 2, 3) of the neutralino in the laboratory 
system by: 

4 = (0,-1,0,0), 4 = (0,0,1,0), 4 = — (g,0,0,-£ x J. (A4) 



Together with p£jm Xi they form an orthonormal set. For the two-body decay x°i ~^ X° n Z 
the decay angle Q\L(j) xv pz) is constrained by sin#?j Iiax = q°/q for q > q°, where q° = 
\2 y m x .,m 2 z ,m' Xn )/2mz is the neutralino momentum if the Z boson is produced at rest. 
In this case there are two solutions 

i^i _ ( m x, +m z -m 2 x Jqcos9 1 ±E Xi ^ \{m 2 Xi) m 2 z , m 2 Xn ) - Aq 2 m\ (1 -cos 2 fli) 

lPzl ~ 2q*(l-cos*0 1 )+2m* i • (A - 5) 

If q° > q, 6\ is not constrained and there is only the physical solution \p%\ left. The 
momenta in the laboratory system are 

Pz = (E z , —\Pz\ sin^! cos 0i, \p% | sin^i sin0x, — |pf | cos0i), (A. 6) 

Pf = (Ef, —\pf\ sin 2 cos 02, \pf\ sin 2 sin 02, —\Pf\ cos0 2 ), (A. 7) 

E '=w = m=m^y (A - 8) 

with 2 l(p Xi ,Pf) and the decay angle D , 2 l(p z ,Pf) given by: 

cos0£) 2 = cos0i cos0 2 + sin0! sin 2 cos(0 2 — 0i). (A. 9) 

The spin vectors t°z (c = 1,2,3) of the Z boson in the laboratory system are chosen by 

^(o.J^f-V 4=(o,ft^L), 4 = -L(|fc|, B *). (A.10) 
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The spin vectors and p z /m z form an orthonormal set. The polarization vectors e Xk,fl for 
helicities = — 1, 0, +1 of the Z boson are defined by: 



8 — ^(^z it 2 z ); e° — t z ; e + — j?(t z + it z ). 



(A.ll) 



B Phase space 



The Lorentz invariant phase space element for the neutralino production (6) and the decay 
chain (7)-(8) can be decomposed into the two-body phase space elements: 



dUps(s,p Xj ,p Xn ,p f ,pf) = 

j—^ dhips(s,p Xi ,p Xj ) ds Xi ^rfLips(s Xi ,p Xn ,p|) ds z dLips(s z ,p f ,p f ), 



± 



(B.1) 



dUps(s,p Xt ,p Xj ) = 
dUps(s Xt ,p Xn ,p^) = 
dUps(s z ,p f ,pf) = 



Q 



1 



sin# d9, 



2(2tt) 2 2|£f gcos^i - E Xi |p| 
dQ 2 , 



dVt\, 



2(2tt) 2 m% 



(B.2) 
(B.3) 
(B.4) 



with s Xi = p Xi , s z = p\ and dQ{ = sin 9i ddi d<pi. We use the narrow width approximation 



for the propagators: / |A(x°)| 2 ds y 



m v - r 



f\A(Z)fds z 



m z r : 



-. The approximation 



is justified for (Y Xi /m Xi ) <C 1, which holds in our case with T Xz < 0(lGeV). 



C Spin matrices 

In the basis (A.ll) the spin matrices J c and the tensor components J cd are 





J 1 



■±- -±- 

V2 V2 




J 2 



J 



22 



J 



23 




J 



J 



33 



J 13 = 




(C.l) 



(C.2) 



.(C.3) 
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Fig. 2a \n\ /GeV Fig. 2b \n\ /GeV 




Figure 2: Contour plots for 2a: cr(e + e —> X1X2)) 2b: BR(%2 — > ZXi)i 2c: <J t = a(e + e — ► 
x B R(X° ^Xi) x BR (^ -»• ^) with BR ( Z -»• ^) = 0.1, 2d: the asymmetry 
*4.£, in the |/z|-M 2 plane for = 0.57T, ^ = 0, taking tan/5 = 10, m = 300 GeV, 
v/i = 800 GeV and (P e -,P e+ ) = (-0.8,0.6). The area A (B) is kinematically forbidden 
by m^o + m^o > y/s {m z + m^o > m%o). In area C of plot 2b: BR(x2 z Xi) = 100%. 
The gray area is excluded by m-± < 104 GeV. 
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matrix elements 



Figure 3: Contour lines of the asym- 
metry At for e+e" -> XxXl\ xl ~> 
Zxi] Z -> ££(£ = e, /i, r), in the ip^-ipMi 
plane for M 2 = 250 GeV and |/z| = 400 
GeV, taking tan/5 = 10, m = 300 
GeV, = 800 GeV and (P e -,P e +) = 
(-0.8,0.6). 
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Figure 4: Dependence on cp^ of the 
vector (Vi) and tensor (Tjj) elements of 
the Z density matrix < p(Z) >, for 
e+e- -> X§ - ^X$, for M 2 = 250 

GeV and = 400 GeV, taking ip^ = 0, 
tan/? = 10, m = 300 GeV, ^/s = 800 
GeV and (P e -,P e +) = (-0.8,0.6). 
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Figure 5: Contour lines of a t = a{e + e- -> xlxV) x BR(x|] -> Z*?) x BR(Z -> 
(5a), and the asymmetry ^ (5b) in the \fi\-M 2 plane for = 0.57T, p^ = 0, taking 
tan/5 = 10, m = 300 GeV, ^ = 800 GeV and (P e -,P e+ ) = (-0.8,0.6). The area A (B) 
is kinematically forbidden by m^o + m^o > ^/s {m z + m^o > m^o). 
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Figure 6: Dependence on ip Ml of the vec- 
tor (Vi) and tensor (Tu) elements of the 
Z density matrix < p(Z) >, for e + e~ — > 
X°X°;X° - ^X?, for M 2 = 250 GeV and 
\fi\ = 400 GeV, taking ip^ = 0, tan/3 = 
10, m = 300 GeV, ^ = 800 GeV and 
(P e -,P e+ ) = (-0.8,0.6). 
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Figure 7: Contour lines of a t 
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Fig. 7b 
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(7a), and the asymmetry At (7b) in the \n\-M 2 plane for ip Ml = 0.57T, ip^ = 0, taking 
tan/3 = 10, m = 300 GeV, y/s = 800 GeV and (P e -,P e +) = (0.8,-0.6). The area A 
(B) is kinematically forbidden by m^o + m^o > y/s (mz + m^o > m^o). The gray area is 
excluded by m-± < 104 GeV. 
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Figure 8: Contour lines of a t 
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Fig. 8b 
Xlxl) x BR(x° 
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(8a), and the asymmetry At (8b) in the |/i|-M 2 plane for ip Ml = 0.57T, ip^ = 0, taking 
tan/3 = 10, m = 300 GeV, = 800 GeV and (P e -,P e +) = (0.8,-0.6). The area A 
(B) is kinematically forbidden by m^o + m^o > \fs (mz + m^o > m^o). The gray area is 
excluded by m~± < 104 GeV. 
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